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We study heat transport in quasi-one-dimensional spin-chain systems by considering the model of 
one-dimensional bosonic spin excitations interacting with three-dimensional phonons and impurities 
in the limit of weak spin-lattice coupling and fast spin excitations. A combined effect of the phonon 
and impurity scatterings yields the following spin-boson thermal conductivity behavior: n a tx T 2 at 
low, k s oc T -1 at intermediate, and k s = const at higher temperatures. Our results agree well with 
the existing experimental data for Sr2Cu03. We predict an unusual dependence on the impurity 
concentration for a number of observables and propose further experiments. 
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Introduction. Recent experiments in one-dimensional 
(ID) quantum magnets have revealed remarkably strong 
thermal transport anomalies associated with the low- 
dimensional spin degrees of freedom. In particular, an 
anisotropic thermal conductivity, comparable in mag- 
nitude to that of metallic systems, was observed in 
quasi- ID spin-^ chain and ladder compounds with a 
large part of the heat current attributed to magnetic 
excitations. P, Only recently, however, a considerable 
progress in the theoretical understanding of the thermal 
transport properties of ID quantum spin systems has 
been made.Q 

In this Letter we consider the heat transport in spin- 
chain systems within the Boltzmann equation framework. 
We develop a microscopic approach to the problem of ID 
spin-boson excitations interacting with the 3D phonon 
environment and impurities. An experimental realization 
of such a system is Sr 2 CuO3.0 

Generally, in order to study transport problems one 
needs to identify a mechanism responsible for the momen- 
tum relaxation. Conventional examples are the Umklapp 
and impurity scatterings. Since the characteristic "band- 
width" for the spin excitations in the real ID spin chains 
materials of interest is very large (J ~ 2000K) in com- 
parison with the experimental temperature range, the 
Umklapp scattering of spin excitations on themselves is 
strongly suppressed (as ~ e~ J / T ). Instead, we propose 
that the leading relaxation mechanism at not too low 
temperatures is the two-stage, bottle-neck process of (i) 
transferring the momentum from spins to phonons, and 
(ii) the subsequent dissipation of the momentum via a 
phonon Umklapp process. The central idea is that the 
stage (i) is the bottle-neck because the spin-lattice cou- 
pling is very weak. Thus, the excess momentum in the 
spin subsystem waits the longest time to get transfered 
to phonons, but once it is transferred it relaxes quickly. 
If this is the case we have to find only the spin-phonon 
relaxation time r sp that now plays the role of transport 
relaxation time for the spin excitations. However, as we 
demonstrate below, the spin-phonon scattering mecha- 
nism alone becomes too weak for the low-lying excitations 



which leads to an infinite spin thermal conductivity k s . 
On the other hand, the impurity scattering is very strong 
in ID at low energies. This scattering mechanism renders 
k s finite for any realistic system and allows us to provide 
a detailed description of the thermal transport. 

Model. The Hamiltonian of a single Heiscnberg anti- 
ferromagnetic spin-i chain is: 

#chain = J ^ Si • Si+l. (1) 

i 

The interchain coupling is very small in Sr 2 Cu0 3 ,l] 
J'/ J ~ 10~ 5 , and is neglected here. After the Jordan- 
Wigner transformation followed by bosonization Eq. l[T)l 
can be expressed as the Hamiltonian of free bosonic field 
<j>, which we refer to as the spin-boson field: 

H = | / dx (n 2 + (d x $) 2 ) = v £ \k\b\b k , (2) 

J k 

where n(a;) is a conjugate field momentum, and velocity 
v is given by v oc J a, a is the chain lattice constant, |5j 
and b, b^ are the spin-boson operators. 

Since the superexchange J is a function of the inter-site 
distance the lattice vibrations can modify it. This will 
couple phonons to the spin-boson Hamiltonian density. 
Such an interaction can be written in terms of the phonon 
operators a, and the spin-boson operators as follows: 

H " = ~% E V,(P,k,k')(al >e blb kl +h.c.) (3) 

^ kk'Pl 

Pwkkl 

V e (P,k,k>) = /a " ,,, (gpi) a Hfe'-fc . ( 4 ) 

where g sp is a spin-phonon coupling constant, the index 
£ enumerates three phonon polarizations, uj-pg = C(\~P\ is 
the phonon dispersion, £ is the polarization vector of a 
phonon, m$ is the mass of the unit cell, N is the number 
of unit cells in the sample, k and k' are the ID momenta 
of spin bosons, and P is the 3D momentum of a phonon. 
Note that only a component of the total momentum along 
the chain is conserved, which is explicitly given by the 5- 
function. Sound velocities ce are much smaller than the 
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spin-boson velocity v for the real spin-chain materials: 
C£ <C v. Eq. J2J describes the processes in which spin 
boson either emits or absorbs a phonon. The projection 
of the phonon polarization vector on the direction of the 
chain (x-axis) (£pe)x is equal to: |PJ/|P| (longitudinal), 
and ^/l — (P||/P) 2 or (transverse), where two answers 
for the transverse phonons correspond to the polarization 
in the plane given by x and P and normal to such a plane, 
respectively. The latter phonons do not couple to spins 
and are not discussed any further. 

The impurity coupling is the final component of the 
model. Impurities lead to a local variation of the su- 
perexchange J, that leads to the scattering of magnetic 
excitations. The most important (relevant in the renor- 
malization group (RG) sense) part of the impurity Hamil- 
tonian for an impurity located at xq can be written as: 

H imp = ^£ cos (2k F x + VMC^(x )^j , (5) 

where the Luttinger-liquid parameter K, = 1/2. 

Phonon mechanism of relaxation. To study the re- 
laxation of ID spin excitations coupled to the bath of 
3D phonons we solve the stationary Boltzmann equation: 
vd x fk = ~~Sk[f], where f k is the spin-boson distribution 
function and S k [f] is the collision integral, given by: 



Sk[f] = f 

Jk> 



w kk ,f k {f kl + 1) - WvkMh + i) 



(6) 



where J,, stands for J dk' /2-k and W kk > is the total prob- 
ability of the spin excitation to be scattered from the 
state k to the state k'. For the processes of scattering 
due to phonons such probabilities are given by: 
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where / p denotes (l/27r) 3 / d 3 P. The probability W k > k = 
W kk i [k <-> k') . We assume that for all relevant temper- 
atures the phonon relaxation time is much shorter than 
that of the spin bosons. Consequently, one can use the 
equilibrium form for the phonon distribution function 
npi = n pe = l/(e Wpf/T - 1). The "elementary" scatter- 
ing probabilities w e kk , p of the spin boson due to emission 
or absorption of the phonon with the momentum P and 
polarization I are determined from Eqs. : 
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where Vo is the elementary cells volume, uj k = v\k\ is the 
spin-boson energy. Note that the total energy and the 
projection of the momentum along the chain are con- 
served, which is enforced by the (5-functions in Eq. 

For a small temperature gradient d x T ^ one can 
introduce a usual ansatz for the distribution function f k 



to linearize the Boltzmann equation: f k = f k +fl, where 
f k is a non-equilibrium correction to the equilibrium Bose 
distribution function /£. Function fl should be small 
and proportional to the temperature gradient. Thus, the 
Boltzmann equation for the spin bosons to the first order 
in d r T can be written as: 



where S k is the linearized collision integral given by S k 
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and 



S/^[k ^ k']. The integral accounts for 
two types of collision events: (i) spin boson with the 
momentum k emits a phonon with the momentum P and 
scatters into the state with the momentum k'] (ii) process 
inverse to (i). Likewise, describes absorption of a 
phonon by the spin boson with the momentum k and the 
corresponding inverse process. 

Once fl is found the spin-boson thermal current den- 
sity J S E can be written as: J% = J, v 2 kfl = —n s d x T, 
where the coefficient k s is the spin-boson thermal con- 
ductivity. Note that the use of Boltzmann equation im- 
plies the semiclassical transport, condition for which can 
be shown to be satisfied for the real ID spin system. 

Kinematic considerations. To solve the Boltzmann 
equation one needs to evaluate integrals in Eq. I|10l) . Part 
of this task can be accomplished without any approxi- 
mations as the integral over the phonon momentum P in 
and can be calculated explicitly. The remain- 
ing integration over k! is restricted to some finite intervals 
due to the conservation of the energy and momentum. 
Further, there is a following simplification. As we 
noted above, without impurities the heat conductivity 
will diverge at small \k\. Thus, the biggest contri- 
bution to k s comes from the low-energy spin bosons: 
v\k\ <§; min(T, Qd) and we need to consider only them. 

Relaxation time. One can write the collision integral 
in Eq. l(TU|) as a sum of two terms: S k = fl/T sp (k) + SS k , 
where the first term has the usual relaxation time form, 
while the second one does not. However, the second term 
can be demonstrated to be small for small spin-boson 
energy SS k /S k oc {vk/Tf. 

Finally, with the help of this approximation, the colli- 
sion integral can be written in the relaxation time form 
Sk ~ f k /iSp(k), where t s ~ 1 (A:) plays the role of the trans- 
port relaxation rate and is given by:Q 

Ak 2 T 3 

Tsvik)- 1 = — — ^ r(0£>/T) , (11) 



where T(z) 



and h{z) = f Q 2 
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the phenomenological constant A is of the order of 
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(yo/mi)(g sp /c) 2 . By evaluating the asymptotic proper- 
ties of T(z) one can easily verify that: 
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for T < e D , 
for T > B D , 



(12) 



where 6_d = ®£>/y/li(oo) w ®d/4 plays the role of a 
crossover temperature. The relaxation rate for the spin 
boson, Eq. (|12f) . is our main result for the spin-phonon 
scattering. We note that the k- and T-dependence of 
the bottle-neck relaxation rate will be the same for other 
low-D spin systems with large J where momentum is con- 
served in the dimensions fewer than D=3. 

It is interesting to note that since the spin bosons are 
fast (c <C v), the energy and momentum conservation 
dictates that the majority of the phonons which interact 
with the spin subsystem must have their momentum al- 
most normal to the chain direction. Indeed, the phonon 
energy is wp = cP = v\\k\ — \k'\\, where k and kl are 
the ID spin-boson momenta, while the momentum con- 
servation along the chain gives |Pm| = |fc — k'\. There- 
fore, for a typical phonon: P ^> |Pii|, except for the case 
of almost elastic backward scattering (k 1 ~ —k) which 
does not contribute substantially to the transport relax- 
ation time. Consequently, we determine £ x i = 0{c/v) 
and £ xt = 0(1). Thus, the most effective spin-phonon 
scattering is due to the transverse phonons. 

With the collision integral in the relaxation time form 
the Boltzmann equation is trivially solved: 

/ i = «!*^ (w ^ R ,_ Z 2M (flLr)j (13) 

and the thermal conductivity is calculated: 
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This expression diverges at small k giving rise to an infi- 
nite thermal conductivity. It occurs because the scatter- 
ing of spin bosons on phonons is not sufficiently strong 
for k — > to ensure the convergence of the low-energy 
contribution to the integral. 

Impurities. In order to remove the divergence of i|14|) 
we need to introduce yet another scattering mechanism 
for the \ow-k spin bosons. In a ID system impurities 
scatter low-energy excitations very effectively. Therefore, 
disorder in the magnetic coupling J will regularize the 
infrared divergence in Eq. I|14l) . Since the spin- impurity 
interaction in Eq. (jSJ) is not a low-order polynomial in <E>, 
the Boltzmann transport theory cannot be directly ap- 
plied to the impurity scattering. Instead, we will evaluate 
the spin boson life-time T lmp using the Green's function 
perturbative expansion in powers of Wimp- 
There are two issues we must clarify before proceeding 
with this approach. First, for the perturbation theory in 
powers of S J; mp to be valid the temperature must be big- 
ger than the Kane-Fisher temperature: T ^ Tkf- For 



our system with K. = 1/2 the Kane-Fisher temperature 
is Tkf = $J\ mp /J- Second, it is important to note that, 
generally, the life-time of an excitation is not equivalent 
to the transport relaxation time. While the only pro- 
cesses which violate the conservation of the momentum 
contribute to the latter, any kind of scattering shortens 
the former. However, since Eq. (JSJ) does not conserve 
momentum and since in any typical scattering event the 
spin-boson momentum changes drastically, there is no 
distinction between these two time scales in our case. 
This justifies the use of our approach. 

The second-order correction to the spin-boson Green's 
function from the impurity scattering is: 



5V[ 2 \r) = 



Bl(r)H inip (r')H inip (T")B kl (0) 



where B\(t) = b\{r) + b_ k (r) and the line over the 
Hamiltonians stands for the disorder averaging. The 
right hand side of this equation is proportional to the 
Matsubara self-energy. After analytical continuation we 
find the retarded self-energy: 



air 2 J\k\T 



(15) 



with n being the dimensionlcss impurity concentration. 
Once the self-energy is found one can obtain the dressed 
Green's function: Dk^ = ^k/{^l — u 2 — 2w/-Sj[! £ w ). By 
finding the poles of this expression the damping of the 
spin bosons is determined as given by: 



T ' 



(16) 



where A 2 oc n5J 2 mp . 7} Unlike in the spin-phonon case 
the impurity relaxation rate is independent of k and thus 
will remove the divergence of the thermal conductivity. 

Thermal conductivity. Finally, having at hands the 
transport relaxation rates for the spin boson due to 
phonon and impurity scattering, Eqs. (jlljl and 116(1 . we 
can calculate the thermal conductivity of the spin chains. 
The total relaxation time can be found according to the 
Matthiessen's rule: r to t 1 = "H 
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expression for the thermal conductivity l|14|) should be 
modified to: 
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4sinh 2 (x/2) x 2 +a(T) 



(17) 



with a(T) = v 5 A 2 /AT 6 T(e D /T). Since the tempera- 
ture is always much smaller than J one can safely replace 
the upper limit in Eq. l|17f) by infinity. 

Eq. H17[) specifies a function with a single maximum at 
T = T m oc (nS J 2 mp / AJ) 1 / 5 , the maximum value of n s is 
K raax _ Ks ^T m ) cx v 6 / AT^. The function k s (T) vanishes 
as T 2 at T T m and saturates at T > T m , approaching 
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FIG. 1: Normalized spin thermal conductivity k 3 vs reduced 
temperature T/T m . Diamonds are the experimental data for 
Sr 2 Cu0 3 , T m = 79.4K and ^ ax = 36.7Wm _1 K _ \ see Ref. 
4. Solid line is the result of this work, Eq. I|170 . for Od ~ 
11.6T m . Arrows mark the saturation value for the solid 
line at T > 8d and the spin-phonon scattering crossover 
scale Qd, respectively. Asymptotic behavior at T < T m (~ 
T 2 ) and at T ^> T m (~ const) are shown by the dashed lines. 
Inset: same data on the log-log scale. 

K oo K K rnax Tm/ Q D . 

f 2.26 T 2 /Tl for T « T m , 
« ^ 1.46 T m /T for T m « T « D , (18) 
* s [5.55T m /6 D forT>0 D . 

At temperatures below T m , k s is controlled by impu- 
rity scattering. At temperature T > T m the number 
of higher energy spin bosons, which scatter mostly on 
phonons, increases. However, the low-fc spin bosons are 
still scattered mostly by impurities. Thus, at T > T m 
both mechanisms contribute to k s . 

Comparison with experiments. Using Eq. (|17|) we com- 
pare our theory with the experimental data for Sr2Cu03. 
It is convenient to measure the temperature in units of 
T m and the thermal conductivity in units of n" lax . This 
guarantees that the location of the maximum coincides 
with the maximum of the experimental curve. Once this 
rescaling is done there is only one parameter to adjust in 
our theory: Q D . By choosing Q D = 11.6T m we fit the 
experimental data for T > T m almost ideally. Shaded 
area at T < T m highlights the region where the phonon 
background subtraction creates a large uncertainty in the 
experimental data.0 With this in mind we conclude that 
Eq. (|17|1 describes very well the whole range of available 
experimental data. 

Our work predicts a remarkable behavior of the ther- 
mal conductivity at higher temperatures: saturation at 
a constant value. A parallel can be drawn with metals, 
where k also saturates at high temperatures. However, 
in our case the excitations are not fermions at Ep, but 
the long- wavelength spin bosons. Because of that, the 
manner in which impurity scattering cuts off such exci- 



tations at low energies remains important even at high 
temperatures. Therefore, in the case of spin chains the 
saturation of result of a non-trivial combination 

of several effects: temperature- and /c-dependence of the 
impurity and phonon relaxation rates and the ID density 
of states of the spin-bosons. We note that the saturated 
regime in k s should persist up to T ~ J above which 
bosonic description of excitations is not adequate. 

Predictions. One of the most straightforward ways to 
verify our theory would be to measure the thermal con- 
ductivity for higher temperatures and check if the satu- 
ration really takes place as we predict. Next, we see that 
the temperature T m at which k s (T) reaches its maximum 
scales as n 1 / 6 . Also, the maximum value of the thermal 
conductivity Kf ax = n s (T m ) scales as n~ 2 / 3 , an unusual 
and rather strong effect. Finally, the saturation value 
has yet different concentration dependence n" 1 / 2 . If 
such a behavior is observed in the materials with the iso- 
tope substitution for either Cu or O ions it will provide 
a strong support to our theory. 

Conclusions. To conclude, we have developed a the- 
ory of anomalous heat transport in the spin-chain sys- 
tems coupled to the 3D phonon environment in the pres- 
ence of weak disorder. We have calculated the thermal 
conductivity as a function of temperature and have ob- 
tained that the low-temperature transport is dominated 
by the impurity scattering while the high temperature 
transport is determined by both the impurity scattering 
and the spin-phonon collisions. Our main results are in a 
very good quantitative agreement with the available ex- 
perimental data. In implementing our approach we have 
also obtained an insight into various microscopic details 
of the problem. This has allowed us to formulate several 
predictions and suggest future experiments. 
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